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The submitted Comment is based on using an isolated quantum dot approach to discuss the situation 
where the coupUng to the leads is considerable (not negligible). This finite lead coupling is the situation 
in most molecular transport junctions. In such situations the population on the molecule is not static and 
dynamical effects (fluctuations) are crucial. While the discussion in Ref. 2, and in the Comment may well 
apply for an isolated molecule, it is simply irrelevant for the very different transport situation analyzed in 
our original contribution. 



o 

(N 

> 

(N 



O 

o 

1 



o 
o 



% 



In their preceding comment on our earlier pa- 
per/ Alexandrov and Bratkovsky (AB) bring 
up the following claims: 

1. A previous work by AB^ that discusses 
bistability in correlated transport through 
a degenerate molecular quantum dot was 
"surprisingly" not cited in our paper.' 

2. The mean field approximation used in 
Ref. 1 is invalid. 

3. The bistability and hysteresis phenomenon 
discussed in Ref. 1 is an artifact of this in- 
valid mean field approximation. 

Indeed, while we have cited the paper on 
bistable tunneling current by AB,^ we indeed 
overlooked the paper 2 by these authors. In 
retrospect we believe (for reasons given below) 
that the paper is only marginally relevant to our 
work, however discussing it at the time would 
have make some issues clearer and could have 
saved the present exchange, so we regret our 
failing to do so. In what follows we address is- 
sues ^ and ^. 

Consider first the mean field approximation. 
Let two coupled subsystems a and b of sys- 
tem s be described by the Hamiltonian Hg = 
Ha + Hh + Vab and let system a be character- 
ized by two states, \la > and \2a >. Assume 
furthermore that an additional process causes 
a rapid interchange between states |la > and 



1 2a >, and that this process by itself would 
bring subsystem a to equilibrium with probabil- 
ity Pi to be in state |la > and P2 to be in state 
\2a >■ How good is an approximation by which 
the dynamics of system b is assumed to be gov- 
erned by the Hamiltonian Hh+ < Vat >a 
where < ... >a denotes an average over the 
state of a, and in particular when can < Vat >a 
be represented by Pi < lalKiblla > +P2 < 

2a\Vab\2a >? 

Describing the motion of subsystem b by 
the Hamiltonian Hh+ < Vab >a (while de- 
riving the dynamics of a from the frozen in- 
stantaneous state of b) is the essence of the 
adiabatic (Born-Oppenheimer) approximation 
used in quantum molecular dynamics, and re- 
lies on the assumption that system a is much 
faster than system b. The correct identifica- 
tion of < Vab >a depends on details of this 
assumption. In particular, in the example con- 
structed above, identifying < Vab >a as given 

by Pi < la\Vab\la > +^2 < 2a\Vab\2a > 

(rather than, e.g. moving system b under the po- 
tential < la|V"afe|la > whcn o is in state \la > 
and under < 2a|T4fe|2a > when a is in state 
\2a >) is valid provided that the dynamics of 
the I la >^^ 1 2a > is fast on the timescale that 
characterizes the dynamics of system b. 

In the present problem system b is the vi- 
bration of frequency wq and system a is the 
molecule, which is characterized by two states 



- occupied and unoccupied. The process that 
interchanges between the two states is the 
molecule-leads electron exchange, whose char- 
acteristic timescale is T^^ = 2itV^p where V 
is the molecule-lead coupling and p - the elec- 
tronic density of states in the lead. From the 
above discussion it is clear that the condition 
for the validity of the mean-field approximation 
as applied in Ref. 1 is F > hujQ. In molecu- 
lar junctions with chemical binding between the 
electrodes and the molecular bridge F is of the 
order 0.1 — 1 eV while typical vibrational fre- 
quencies span the range 0.001 — 0.3 eV. There is 
thus a wide regime in which this approximation 
is valid. 

In the preceding comment, AB illustrate the 
failure of the mean field approximation for a 
model in which F = 0. Clearly this illustra- 
tion is irrelevant for the problem at hand. They 
then proceed to discuss, under what they call 
"the correct procedure", a different model (a re- 
cast of Ref. 2) in which bistability stems from 
electron-electron interaction in a d-fold degen- 
erate molecule. Again, analysis of this model 
is done by neglecting the molecule-leads cou- 
pling and, furthermore, disregarding the effect 
of the electron-vibration interaction on the state 
of the vibration - two principal ingredients of 
the model of Ref. 1. While we have no ob- 
jections to this model and its consequences, it 
can hardly be regarded as a "correct procedure" 
with respect to our model, or, indeed, to the 
molecular junction problem that we analyzed. 
We note that similar models, with similar con- 
clusions about the inherent bistability in the 
system (though with different views about its 
consequence - see below) have been recently 
discussed by several authors.'*"' 

Next consider the consequences of the 
bistable nature of the nuclear potential surface 
obtained in the mean field approximation. Ob- 
viously, as emphasized by Mitra et al,^ this 
bistability does not imply any phase transition 
property. The observation of hysteresis or, in- 
stead, switching behavior (telegraphic noise) 
depends on the timescale of changing the ex- 
ternal control parameter (e.g. the bias or the 



gate potential) relative to the rates of transitions 
between locally stable states.' Both phenom- 
ena are observed in molecular junctions, and 
the question whether they are accounted for by 
the present mechanism should be settled by de- 
termining these rates. This issue has not been 
raised or discussed in the preceding comment 
byAB. 

We end our reply with several comments: 

1. For an isolated molecule (MQD in the 
AB language) the molecular level is well de- 
fined, and is shifted from cq to cq — NP /ujq 
by the electron-vibration interaction (reorgani- 
zation energy). Coupling to the leads gives rise 
to electronic level population that fluctuates on 
the characteristic timescale F^^. The observa- 
tion that for F > Tiujq the electronic energy 
shift does depend on the average population of 
molecular level, was first made in Refs. 8, 9. 
The isolated MQD consideration of 2 misses 
the physics of this situation, which is physically 
relevant to molecular transport junctions. 

2. The AB Comment refers to the approxima- 
tion v? = riQ-fi as "a spurious self-interaction 
of a single polaron with itself". In fact, this 
term represents the interaction of a tunneling 
electron at a particular time with the polariza- 
tion charge established in response to electrons 
(many tunneling events) which traversed the 
junction at earlier times. This "self-interaction" 
is legitimate in the (physically relevant) case 
where coupling to the leads (disregarded in the 
AB treatment) is taken into account. 

3. The treatment of Ref. 7 of the same one- 
level model, which goes beyond the mean-field 
approximation, also leads to multistability in 
agreement with our treatment. The statement 
in Ref. 2 that "the Born-Oppenheimer approx- 
imation does not apply to nondegenerate level, 
since there are no fast (compared to the charac- 
teristic phonon time 1/ujq) electron transitions 
within the dot" holds for the isolated MQD. 
In molecular wires the molecule-metal electron 
exchange with timescale F~^ provides this fast 
process. 

4. Contrary to statements made in the AB Com- 
ment, taking into account the coupling with the 



leads (after the Lang-Firsov transformation ap- 
plied to the Hamiltonian as is done in Ref. 2) 
does provide nonlinearity, since the effective 
couplings do depend on the electron population 
on the bridge. Confusion here is due to the ap- 
proximation made in 2 which neglects the effect 
of the electron- vibration coupling on the vibra- 
tion state. For detailed discussion see 10. 

In conclusion, while the discussion in Ref. 2, 
and in the AB Comment apply for an isolated 
molecule, it is simply irrelevant for the very dif- 
ferent transport situation analyzed in our origi- 
nal contribution. 
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